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Abstract 

A realistic model is suggested based on the quasiparticle Lagrange version of the self- 
consistent Finite Fermi Systems theory supplemented with the microscopically calculated sur¬ 
face parameters of the Landau-Migdal interaction amplitude. The latter are expressed in terms 
of the off-shell T-matrix of free iViV-scattering and show a strong dependence on the chemical 
potential of a nucleus under consideration in the drip line vicinity. This effect could result in 
shifting the neutron drip line position to very large values of the neutron excess. 
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Up to now, all predictions on the location of the nuclear drip line for heavy nuclei are based 
on phenomenological approaches. All of them operate with parameters adjusted to the prop¬ 
erties of stable nuclei. This point seems to be questionable for nuclei far from the /^-stability 
valley. In view of commonly recognized importance of pairing for the drip line, the main efforts 
in this field were focused to study the superfluidity effects at small values of the chemical poten¬ 
tial fi BBi- Such problems were investigated as a correct account of continuum, comparison 
of exact Bogolyubov solutions versus those within the BCS approximation, the surface pairing 
versus the volume one, and so on. In this letter we concentrate on the examination of the 
neutron average potential well at small values of fin- We give arguments in favour of significant 
variation, in the vicinity of the drip line, of some parameters of the effective A^A^-interaction 
which generates this potential well. It could result in shifting the position of the drip line into 
the region of large (or, maybe, very large) values of the neutron excess. 

The reason for such an effect can be readily shown in terms of the simplest version of the 
self-consistent (SC) Finite Fermi Systems (FFS) theory ii which is based on the simplified 
version of the self-consistency relation of Ref. 1^: 



where F is the Landau-Migdal (LM) amplitude. Obvious isotopic indices in Eq. (^ are omitted. 
Let us for a while limit ourselves to the zero-range components of F, which in the standard 
FFS theory notation has the form: 


F(r, r) = Co [/o(r) + /o(r)rir2 + (^o(r) + ^Ur)UT-2) (Ti(T2\ (5(r - r') 


( 2 ) 


where the normalization factor Co = {dn/dep)~^ is the inverse density of states at the Fermi 
surface. 

In the FFS theory, a strong dependence of the scalar-isoscalar amplitude /o on the observa¬ 
tion point r was introduced. In fact, in Ref. @] a simple interpolation form of such a dependence 
was suggested: 

,P+(r) 


fo{r) = r + {r-ry- 


Po 


( 3 ) 


Here p+(r) = Pn(r) -|- Pp{r) is the nuclear density in the point r, while po = p+(r=0). The 
subscript ”0” for the zero Landau harmonics is for brevity omitted in the r.h.s. of Eq. (^ and 
below. 

It is worth to mention that the density dependence of the phenomenological Skyrme force 
0 agrees with the ansatz of Eq. (^). There exist also alternative versions of the interpolation 
formula for /o(r) by means of the function {p/po)°‘,a ^ 1, in the r.h.s. of Eq. (^) or a more 
complicated density dependence p|. 
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. But all of them are characterized by a strong 
difference between the dimensionless parameters and /™: ~ —3, whereas is close to 

zero. [] 

It should be noted that the density dependent scalar-isoscalar amplitude 


/'(r) = /'--F (/'--/'-) 


P+(r) 

Po 


( 4 ) 


^The inequality /“ > — 0.5 should be fulfilled to avoid the Pomeranchuck instability 
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was used in a new version of the energy functional method by S.Fayans et ah p, based on a 
detailed analysis of long isotopic chains. The difference between the parameters is also 

signihcant, though not so dramatic as in the isoscalar case. In principle, relations similar to 
Eq. d^) can be written also for the spin-dependent terms of Eq. (|^), but up to now no evidence 
of a noticeable difference between the internal and external values of these amplitudes was 
found. Therefore the equalities g'^^ = g'^^ were imposed in the FES theory. 

In Ref. it was found that the external values of the LM amplitudes can be calculated 


in terms of the off-shell T-matrix of free A^iV-scattering taken at a negative energy E = 2/x, 
where /r is the chemical potential of the nucleus under consideration. Stable nuclei with /i„ = 
/ip ~ —8 MeV were considered, and a reasonable agreement with the phenomenological values 
for the surface parameters of the LM amplitude was obtained. 0 It conhrmed the relevance of 
the asymptotic relation F —> T(2/i) for the description of the properties of stable nuclei. 

For the spin independent amplitudes under consideration the explicit form of these relations 
is as follows: 


/o“ = ^ [‘o(B = 2 a) + HE = 2a)| , 
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/r= [‘o(B =2a) 


3ti(E = 2a)] , 


( 8 ) 

( 6 ) 


where to, ti are the dimensionless values of the off-shell T-matrix with the spin value S' = 0 and 
S' = 1 respectively, taken at the zero value of all nucleon momenta. In stable nuclei the isospin 
symmetry works well and one has: 


= fr + /r = \h(E = 2 a), (7) 

fZ = /o“ - /r = 5 l*o(-E = 2 a) + 3(.(B = 2a)] , (8) 

In this paper these relations are applied to nuclei close to the drip line in which the neutron 
and proton chemical potentials deviate signihcantly from each other. As far as both amplitudes 
to and ti at small energy E depend on E rather sharply, the isospin symmetry is destroyed. 
As a result, one deals with the situation when /“ 7 ^ Both of these amplitudes can be 
found from Eq. (J^), but at T = 2/i„ for the hrst one and at T = 2/ip, for the second one. The 
amplitude /®p can be obtained from Eq. (||) at T = /i„ -|- /ip. In the neutron drip line vicinity, 
only the neutron chemical potential is close to zero, whereas, on the contrary, |/ip| increases 
with approaching the boundary. Therefore only the amplitude /“ changes signihcantly (see 
Fig. 1). 

The SC relation, Eq. (|l|), with the amplitude F given by Eqs. (@), (Q), (H) can be readily 
integrated over for a spherical nucleus yielding the following relation for the neutron potential: 


UJr) = Co 


fTn Pn(r) + (/” - /S)^ (pAt) + pH)) 
+fZ pH) + U'Z - /S)^ (/’+(’■) + PvH) 


(9) 


^They depend a little on the type of the density dependence. The best agreement was achieved with the 
spin independent parameters of Ref. and the spin dependent parameters of Ref. 0 
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A similar relation is obtained for the proton potential Up{r) replacing the indices ”n” and ”p”. 

To estimate the effects of the energy (or p-) dependence of the interaction amplitude /“, 
let us hnd the value of Un{0). To simplify the expression, we use the approximations Pn(0) = 
{N/A)pq and Pp(0) = {Z/A)pQ. Then we get: 


Un{0) 


-Copo 


/, 


ex I £ 1 ] 

nn ' Jn 



+ fi 


ex- f 

J^2 


np 



( 10 ) 


Let us add a small number of neutrons to a heavy nucleus. Then, neglecting for a while the 
pairing effects, we obtain an approximate relation 

6pn = SU^iO) ( 11 ) 

for the change of the neutron chemical potential. | The variation of the expression ( 0 ) is a 
sum of two terms, 

6Un{0) = SUi + 6U2 , ( 12 ) 

having different origin. The hrst term results from variation of N and A while the second one, 

m = IcoPo^^/S . (13) 

is due to the //-dependence of the amplitude Neglecting this //-dependence one gets the 
accustomed variation of the chemical potential 

6^1 = 6U,. (14) 


But taking of the second term into account yields a dramatic deviation from this traditional 
result when nuclei with a big neutron excess are considered. 

Indeed, at small values of //„, the amplitude /“, Eq. ( 0 ). taken at E = 2//„ is singular: 


/ ex _ 

nn 


a 


a 


y/E yfPp^r 


(15) 


Here we neglected the small, compared to //„, value of the virtual level energy of the T-matrix 
in the singlet channel. The variation of Eq. ([T5|) yields: 





(16) 


Upon substituting Eqs. 


([I^) , and (pul) into Eq. (|^) one hnds 




1 + ho/ {2pn) 


(17) 


where a short notation Vq = for the hrst term of Eq. (|I0|) is introduced. 

^This equation contains also an additional term provided the neutrons occupy a new j-lGvel. All the consid¬ 
eration presented below remains valid in this case, too. 
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It can be easily seen that the denominator of the relation (|^) is noticeably greater than 
nnit for nnclei with small value of Let us, for example, calculate this quantity for two 
isotopes of tin in the vicinity of the old drip line {A^ax = 176 which is the common value to 
all calculations for the tin isotopes as far as we know). We use the values of /i„ which will be 
found below and the corresponding values of the amplitude taken from Fig.l. We take also 
the standard values of the normalization parameters: Cq = 300 MeV-fm^ and po = 0.16 fm“^. 
First, let us consider the ^®°Sn nucleus. In this case, we have = —1.4, Vq = —15 MeV, 
and pn = —3.4 MeV. The substitution of these values into Eq. (p!?]) yields 6pn = ^hn/3-2. The 
analogous calculation for the ^°°Sn isotope {pn = —2.0 MeV, = —1.66, Vq = —22.4 MeV) 
results in 6pn = 6p^/6.5. We see that deviations from the traditional scheme are really large and 
are growing as soon as with the value of \pn\ becomes less and less. This explains qualitatively 
why nuclei (e.g. ^°°Sn ) which are strongly unbound in traditional calculations could exist in 
our approach. 

Going to the actual calculations incorporating the effect discussed above, we start from an 
advanced version of the SC FFS theory, the so-called quasiparticle Lagrange method (QLM). Q 
It was devised in Ref. for magic nuclei and extended in Refs. [0, to include superfluidity. 
All necessary modihcations in the scheme can be explained for a simpler, nonsuperfluid case. We 
present here only a brief sketch of the QLM referring to Refs. |^, for details. The approach 
utilizes the Lagrange formalism which is the most convenient when the energy dependence 
effects are considered explicitly. The effective Lagrangian is constructed in such a way that its 
variation results in the Dyson equation. 


(£-£?-E,)G, = 1, 


(18) 


for the quasiparticle Green function Gq with the quasiparticle mass operator 

Sg(r,k,£) = So(r) + -^ k Si(r) k + ^ S 2 (r), (19) 

fcp 

where the normalization parameters are: fcp = 7r/(mC'o), £p = (/cp)^/(2m). 

The hrst two terms of are common to the HF theory with effective velocity dependent 
forces (e.g., Skyrme forces). The third term, S 2 (r), is due to energy dependence. It determines 
the coordinate dependent Z-factor: 


Z{t) = (l - Ej(r)/£») 


-1 


( 20 ) 


The solution of the Eq. (p!8[) can be expressed in terms of the eigenfunctions 4/^ of the 
corresponding homogeneous equation as: 


G'g(ri,r2,£) = ^ 


T^(ri)TA(r2) 


e-ex + i6sgn(ex - p) ’ 
where 5 is small and positive. The functions 4 /a are orthonormalyzed with the weight: 


J dr^l(i)Z '(r)'I'y(r) = iiy. 


( 21 ) 


( 22 ) 


^The term ’’quasiparticle” is used here in Landau’s (not Bogolyubov’s) sense. 
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The quasiparticle density z^o(r) associated with the TA-functions is: 

i^o(r) = ^uaITaWI^ (23) 

A 

where ua = (0, 1) are the quasiparticle occupation numbers. It differs from the usual density 
p(r) normalized to the total particle number by the Z-factor: 

z/o(r) = Z(r)p(r). (24) 


There are two additional densities introduced in Ref. [Q, the quasiparticle kinetic energy 
density 


1 


i/i r = 




^nA|VTA(r)|2, 


and the total quasiparticle energy density 




(25) 


(26) 


'F A 


The density i^i(r) is analogous to the quantity r(r) of the HF theory [0, whereas the density 
1^2 (r) is a new ingredient of the QLM, which does not appear in the HF approach. 

In notation of Ref. the density of the interaction Lagrangian L' is 


C'(t) = -Cn 


+ ^^i^o+{T)uon{T)uop{r) + Aoii/o+(r)z/i+(r) 


+Aoii/o-(r)z/i_(r) + Ao2i^o+(r)z/2+(r) 


AooFq 


(Vz/o+(r))^ 


+ ^1, 


(27) 


where the normalization density is Pq = 2 (/cp)^/( 37 r^), and Ci includes the Coulomb and the 
spin-dependent terms (mainly, the spin-orbit one). It is worth to mention that the isotopic 
structure of the ’’triple” term in Eq. (0 (proportional to 7 ) is similar to that of the Skyrme 
Hamiltonian of Ref. 0 . 

To extract from the Lagrangian L' the LM amplitude F, which, in accordance with the 
Landau prescription, is the second derivative of the ground state energy, one needs to hnd the 
corresponding Hamiltonian expressed in terms of the usual densities p(r) and r(r). Then it 
is necessary to calculate the second variational derivative with respect to p. This leads to a 
rather cumbersome expression which can be found in Ref. |^. However, the relation of the 
constants Aqo, Aqq with the external LM amplitudes in Eqs. (|^), (^) can be understood 

without using any explicit formula. Indeed, in the asymptotic region outside the nucleus the 
Z-factor is tending to unit and the densities r'o(F) and p(r) coincide. Therefore the variation of 
V with respect to p can be replaced by the variation with respect to z/q. The amplitudes under 
consideration originate from the first three terms of Eq. (^). But the second variation of the 
third term vanishes at large r and one ends up with the following identities: 


AOO — Aqq — / 


Of course, they follow also from the explicit relation for F of Ref. 
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( 28 ) 
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Dealing with superfluid nuclei, we utilize the modiflcation of QLM developed for this case 
in Refs. |1l3| , |1^]. Although the method itself is, in principle, rather general, the practical 

is quite simple. The main approximations are as follows. First, the A- 

£mm = -(20^-25) MeV, = 


'mm ^ ^ ^max5 


scheme of Refs. [0, | 

representation with a limited A-basis is used 
5 MeV), with the discretization of the continuum. Second, the (5-form density independent 
(’’volume”) pairing interaction is considered with the strength depending on the basis, in accor¬ 
dance with the prescription of Ref. |^: F^ = (Fq ln“^(C'p/^), where ^ = \J{jJi — £mm)(e^max — /w)- 
At last, the diagonal approximation for the gap A is used: Aaa' = Thus, the method 

of describing the pairing effects in Refs. [|T^ 113 in the main points coincides with the BCS 


approximation. Although such a scheme possesses some well-known deficiencies for nuclei near 
the drip line and there exist much more advanced approaches for this case (see, e.g.. Refs. @,0), 
we conserve here all details of the pairing scheme |^, (including the values of the pairing 
parameters) in order to separate the effect of the energy dependence of the mean held more 
clearly. 

Our goal is to single out effects of the energy dependence of the external values of the LM 
amplitude. For convenience, let us rewrite the first three terms of Eq. 0 ) substituting the 
quantities /®^, for Aqo, Aqq: 




= -On 


f 


f 


/ex 
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— i^o+{r) + — z/o_(r) + ^ i^o+(r)', i^on(r) z/op(r) 


(29) 


If we consider as phenomenological parameters, there is no difference between this 


expression and the initial one, just the physical meaning of the constants Aqo, A( 


00 


IS more 


transparent. The next step is the calculation of these parameters from Eqs. (i, i)- They 
are = — 2.6 and = 1.56, Ref. |^, instead of Aqo = —3.25 and Aqq = 2.4 of Ref. [§ 
^ Now Eq. (^^ contains only one adjustable parameter 7 instead of the three those of the 
corresponding part of Eq. (EZD- We chose it in such a way to better reproduce the single¬ 


particle energy spectrum of the nucleus ^^^Sn calculated in Ref. |]^ where it was used for fitting 
the parameters of the calculation scheme and was found in a reasonable agreement with the 
experimental one. It should be noted that we do not analyze here the total binding energies and 
the density distributions but concentrate on the single-particle spectra because they are mainly 
responsible for the position of the drip line. We found that the spectrum obtained with the ab 
initio constants and with the value of 7 = 1.6 (instead of 7 = 3.2 in Refs. |^, |^) is in 

a reasonable agreement with that of Ref. [Q (see Table 1), where E\ = (e^ — A^)^ + A^. 

The generalization of Eq. (^) to nuclei with a large neutron excess, where the isotopic 
symmetry is violated, is quite obvious: 


Anfr) = C-n 


^/nn(^ = 2/in) = 2/ip) ^Ipir) 

27 


+/np(^ = /in + /ip) Zi0n(r)zi0p(r) + «io+(r) iion(r) iiop(r) 

'jPo 


(30) 


where the external interaction LM amplitudes Znn; fppi fnp energy dependent and should be 
calculated for the nucleus under consideration in the same way as in Eq. (^. Again the energy 

®These values of /®^, /'®^ are closer to those of Ref. |^. 
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dependence of the first of them only is essential in the case under consideration. In the two 
next amplitudes, this dependence is retained just for presenting a more general form which is 
essential, e.g., for nuclei near the proton drip line. 

We made a series of the self-consistent calculations for the chain of the tin isotopes, hrst, 
with the phenomenological Lagrangian, Eq. (1^) , ^ and, second, with the semi-microscopic 
Eq. (|30|) . Results are displayed in Fig. 2, together with the predictions of Dobaczewski 

The latter two approaches, just as the one of Ref. |T^ 


one 
et ah 


and Fayans et ah 


use 


the phenomenological energy-independent forces Q The results of the three phenomenological 
calculations are quite close to each other. Qualitatively, they could be considered as one 
’’phenomenological” curve which is in very good agreement with available experimental data. 
At small and intermediate values of the asymmetry parameter y = {N — Z)/A the results of our 
semi-microscopic calculation deviate from the phenomenological curve (and, consequently, from 
the experimental data), but not signihcantly. The reason of this dehciency is quite simple. The 


matter is that the one-parameter Eq. (pOf) does not permit to obtain simultaneously reasonable 
values of the two ’’inner” amplitudes, /“ and /'™. It is well known that the Skyrme prescription 
of the isotopic structure of the ’’triple” term used in Eq. (|30D is not obligatory for the effective 
force. The simplest two-parameter generalization is: 


3p§ 


t^o+(r) z/on(r) 


6po L 


-fvl+ir) + ivQ+{v)pl_{r) 


(31) 


which reduces to Eq. (pOD at 7 ' = — 7 . The ansatz of Eq. ( |3TD seems to be more adequate to 
ideas underlying the approach suggested. Such a generalization, with a new adjustment of the 
parameters, will be carried out in a separate work. Preliminary estimates show that in this case 
the ” semi-microscopic” results for stable nuclei become much closer to the phenomenological 
ones. 

When we approach to the drip line, deviations grow, the semi-microscopic curve being sig¬ 
nihcantly higher than the phenomenological one. It is remarkable that our calculation predicts 
the existence of nuclei beyond the commonly admitted end of the tin chain, = 176. It is 
worth to mention that this value of Amax corresponds to the value oiy = 0.43 which exceeds the 

m 


critical value j/o = 0.37 predicted for the asymmetric nuclear matter 


The energy depen¬ 


dence effects under consideration make this difference between hnite nuclei and nuclear matter 
more pronounced. We interrupted our calculations at A = 208 just because we are conscious 
that the dehciencies of the approach become more serious with the value of |/in| becoming less. 
One of them is too schematic consideration of pairing. Then, corrections to the standard FES 
theory could appear when the energy dependence of the LM amplitude is signihcant. 

It is worth to note that the two recipes of taking into account the energy dependence 
effects, Eq. (|j) and Eq. (^0]), are not identical. Indeed, in the hrst one the amplitudes /“,/'“ 
are completely energy independent. As to the second recipe, the parameter 7 is supposed 
to be energy independent that leads to some energy dependence of the quantities /™, /'“ 
evaluated via the procedure described above. In the framework of a pure phenomenology, it is 
difficult to choose between the two ansatz. In our regular calculations, we use the second choice 


®They just repeat the calculations of Ref. Q 

^The trivial, linear energy dependence, taken into account in Ref. Q 
via the density V 2 (j), is not important for the effect under consideration 


which is incorporated into Eq. (UT\) 










just because we have a well developed calculation scheme for this case. A more consistent 
microscopic theory for solving this problem is necessary. Such an approach based on the 
Brueckner theory for hnite nuclei is now in progress |^, |^, but it is yet far from to be 
completed. A semimicroscopic model for the scalar-isoscalar LM amplitude /(r) suggested in 
Ref. on the base of an approximate calculation of the G-matrix for a slab of nuclear matter 


17| shows some energy dependence of the parameter /™, but it is rather smooth. Evidently, 


the truth is somewhere between the two ansatz discussed above. What is important for our 
preliminary analysis, that both of them predict qualitatively the same. Namely, the energy 
dependence of the external LM amplitude /“ makes the neutron potential well deeper when 
the absolute value of fin becomes less. This effect could help nuclei with a huge neutron excess 
to survive. 

The calculation scheme proposed in this letter possesses some evident dehciencies partially 
discussed above. Improvement of them could make the effect less pronounced. That is the 
reason why the title of the paper is formulated as a question. However, the effect itself is based 
on a transparent physical phenomenon. Therefore we expect that it should remain, at least 
qualitatively, after all necessary corrections to the approximate calculation presented in this 
letter. 

This research was partially supported by Grants No. 00-15-96590 and No. 00-02-17319 from 
the Russian Foundation for Basic Research. We thank S.V. Tolokonnikov for useful discussions 
and for supplying us with the numerical data used in Ref. P]. 
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Table 1 


Single-particle spectrum of the nucleus. 


A 

Ex, Mev 

Ex, Mev 


11 

this work 


m 

this work 

exp. 

2pi/2 

-16.86 

-17.94 

-16.96 

-18.05 


l5'9/2 

-15.15 

-15.93 

-15.29 

-16.11 


2^5/2 

-10.19 

-10.82 

-10.42 

-11.05 


3Si/2 

-9.53 

-8.82 

— 

9.97 

-9.18 

-8.64 

l5'7/2 

-8.34 

-8.32 

— 

8.81 

-9.19 

-9.63 

2(i3/2 

-8.08 

-8.05 

— 

8.68 

-8.66 

-8.52 

lhii/2 

-7.19 

-6.87 

— 

5.83 

-5.08 

-5.73 

2 / 7/2 

-2.34 

-2.29 

— 

2.24 

-2.14 


3p3/2 

-0.97 

-0.96 

— 

0.94 

-0.91 


3pi/2 

-0.42 

-0.26 

— 

0.39 

-0.23 
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Figure captions 


1. External LM amplitudes f^[E), and f^{E) (the latter is divided by 10) taken 

with negative sign. The solid lines show these amplitudes within the energy limits which 
are typical for tin isotopes. 

2. A half of two-neutron separation energy S'2n/2 calculated for even-even tin isotopes. Re¬ 
sults of the semi-microscopic calculations of this work (solid line) are displayed together 
with predictions of Ref. |T^ (dashed line), Ref. @] (triangles) and Ref. 0 (squares). 
Experimental values of S'2n/2 are shown by circles. 
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